Evidence for differentiation in the iron-helicoidal-chain in GdFe 3 (B0 3 ) 4 
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We report on a single-crystal X-ray structure study of GdFe3(BOs)4 at room temperature and at 
T=90 K. At room temperature GdFe3(BOs)4 crystallizes in a trigonal space group R32 (No. 155), 
the same as found for other members of iron-borate family RFe3(B03)4. At 90 K the structure 
of GdFe3(BOs)4 transforms to the space group P3i21 (No. 152). The low-temperature structure 
determination gives new insight into the weakly first-order structural phase transition at 156 K and 
into the related Raman phonon anomalies. The presence of two inequivalent iron chains in the 
low temperature structure provides a new perspective on the interpretation of the low-temperature 
magnetic properties. 

I. INTRODUCTION 



£f~) [ The family of borates RM3 (603)4 with a rare earth (RE) or yttrium as R and Al, Ga, Fe, or Sc as M crystallize 
in the huntite, CaMg3(C03)4, structure type with space group i?33i*2iii. The interest in this family of crystals arises 
both from a fundamental point of view and from already realized and proposed applications. Crystals of YAI3 (603)4 
and GdAl3(B03)4 doped with Nd were widely studied during recent years and have been used in optical devices, such 
D ' as self- frequency doubling and self-frequency summing lasers (see, e.g., Refji and references therein). Concentrated 
5_i , NdAi3 (603)4 crystals are efficient media for minilasers^. 

Apart from the interesting optical properties arising chiefly from the lack of inversion symmetry, the 'sub-family' of 
j * I RFe3(B03)4 also attracts considerable attention due to their structure peculiarities. The presence of magnetic order 
at temperatures less then 37 K was attributed to magnetic Fe-Fe or Fe-O-Fe interactions inside quasi one-dimensional 
(ID) iron chains 3 ^. Recent works&£, focusing on low temperature magnetism in GdFe3(B03) 4 , revealed two magnetic 
phase transitions. The second-order magnetic ordering phase transition at Tni =37 K (antiferromagnetic ordering of 
Fe atoms) is followed by a first-order spin-reorientational phase transition at Tn2 =10 K. Additionally specific heat 
and Raman measurements on single crystals^ of GdFe3 (603)4 revealed a weakly first-order structural phase transition 
q ■ at T s =156 K. This structural phase transition is observed in almost all members of the RE ferro-borates familyS&S. 
Recently, Hinatsu et al£ have shown that DyFe3 (603)4 undergoes such structural phase transition at 340 K, by 
■ measuring the temperature dependence of the lattice parameters on powder samples. Here too, a peak in the specific 



heat was observed at this temperature. Similar peaks in specific heat vs temperature dependence were found for 
other heavy RE (R= Eu-Ho, Gd, Tb) ferroborates. These peaks were ascribed to structural phase transitions. The 
£SJ ' transition temperatures T s were found to depend linearly on the ionic radius of the RE. To date, there are no single 
. crystal data available nor has the low-temperature (LT) space group been determined. 

The high temperature (HT) i?32 structure of powder RFe3 (603)4 compounds was first determined in Ref. 1 for R 
= La, Nd, Sm-Ho, and Y. X-ray experiments on single crystals with R=Ndo.5Bio.g^ and R=La, Nd, and Yo.56io.g2 
confirmed this structure. Moreover, Raman measurements for LaFe3 (603)4^, NdFe3 (603)4^^ and GdFe3 (603)4^ 
are in good agreement with group theoretical analysis, based on the i?32 structure. 

Summarizing, whereas the HT crystal structure of most of the ferroborates is known, the LT space group, lattice 
parameters, and atomic positions are unknown. In this work we report on a single-crystal X-ray diffraction struc- 
ture study of gadolinium-iron borate at room temperature (RT) and at 90 K . It is found that the LT-structure 
has the P3i21 symmetry. Two nonequivalent spin chains exist in the LT-structure which gives new insight in the 
low-temperature magnetic properties. The P3i21 symmetry of the LT-phase also confirms the interpretation of the re- 
ported IR absorption by crystal-field transitions^*^ and leads to a better understanding of some peculiarities observed 
in Raman spectrai*^. 



II. EXPERIMENT AND RESULT OF STRUCTURE DETERMINATION 

Crystals of GdFe3 (603)4 were grown using a K2M03O10 - based flux, as described in RcfA 6ig transparent single 
crystals of gadolinium ferroborates were light green in color and had a good optical quality. A block-shaped crystal 
('broken- fragment') with the dimensions of 0.22 x 0.15 x 0.11 mm 3 was mounted on top of a glass fiber and aligned 
on a eruke»i£ SMART APEX CCD diffractometer. The crystal was cooled to 90(1) K using a 8ruker KRYOFLEX. 
Intensity measurements were performed using graphite monochromated Mo-Ka radiation. Generator settings were 



50 KV/ 40 mA. SMART— was used for preliminary determination of the unit cell constants and data collection 
control. The intensities of reflections of a hemisphere were collected by a combination of 6 sets of exposures (frames). 
Each set had a different angle for the crystal and each exposure covered a range of 0.3° in u>. A total of 3600 frames 
were collected with an exposure time of 10.0 seconds per frame. The overall data collection time was 16 hours. Data 
integration and global cell refinement was performed with the program SAINT. The final unit cell was obtained from 
the xyz centroids of 4767 and 5439 reflections after integration for RT and 90 K, respectively. Intensity data were 
corrected for Lorentz and polarization effects, scale variation, for decay and absorption (an analytical absorption 
correction was applied), and reduced to F%. The program suite SHELXTL was used for space group determination 
(XPREP). 




9 n _Fe ~-~-4>-> -B0 3 
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FIG. 1: The structure of GdFe3(BOs)4 in two different projection. The left panels (a) and (b) show the RT structure, and the 
right panels (c) and (d) show the structure at 90 K. The Fe atoms are arranged in chiral chains parallel to the c-axis. Different 
chains are separated by GdOe prisms and BO3 groups. The unit cell outline for R32 is shifted by (1/3, 1/3, 0) for comparison 
with the LT P3i21 structure. 

At both temperatures of RT and 90 K the unit cell was identified as trigonal: reduced cell calculations did not 
indicate any higher metric lattice symmetrjii^. Space groups i?32 and P3i21 were derived for RT and for 90 K, 
respectively, from the systematic extinctions and discriminated from other candidate space groups, which comply 
with the same extinctions conditions, during the structure determination process. Examination of the final atomic 
coordinates of the structure did not yield extra crystallographic or metric symmetry elements^i^. 

The polarity of the structure of the crystal actually chosen was determined by Flack's x-refinemenliiSiiLiSiiS; 
refinement resulted in a x value of 0.50(1), so ultimately an inversion twin is used in the refinement. 

Crystal data and numerical details on data collection and refinement are given in Table [I] Final fractional atomic 
coordinates, equivalent displacement parameters and anisotropic displacement parameters are given in Table [H] for 
data at RT and in Table UTTI for the 90 K data. 



III. DISCUSSION OF THE STRUCTURE 



After a brief introduction describing the main features of GdFe3(B03)4, we will focus on the difference between the 
RT and the 90 K (LT) structure. First, we will try to understand the consequences that the structural changes have 
on the interpretation of Raman spectra anomalies at the weak first-order phase transition reported in the literature!. 
Second, we will give a new perspective for the interpretation of low temperature magnetic data&£. 

The room temperature structure of GdFe3(B03)4 belongs to the R32 space group (Fig. ^ and^)). Our measure- 
ment confirms the structure reported previously for powder samples^ and for several single crystals^ from the family 
RFe3(BOs)4. The structure consists of alternating layers (parallel to afc-plane) of Fe-Gd and BO3 groups (see Fig. 
lb). The main features of this structure are already described in literature^: the BO3 groups are arranged in layers 



TABLE I: Crystallographic parameters 



a-Crystal data and details of the structure determination. 

Room Temperature = 293(1) K 



Moiety Formula 
Formula Weight (g.' 
Crystal System 
Space group. 
a (A) 

c (A) 
V (A 3 ) 

(-) range unit cclhmi 
deg; reflections 



Fori 



Space group_Z 

Z' ( = Formula_Z/SpaceGroup_Z) 

F(000) electrons 
fi( M o K a) (cm ^ ) 

b-Data collection. 

\(MoKa), A 
Monochromator 
Measurements device type 

Detector area resolution (pixels/mm) 
Measurement method 
8 range; min. max., deg 
Index ranges 

Min.- Max. absorption transmission factor 
X-ray exposure time, h 
Total data 
Unique data 

Data with criterion: (F Q > 4.0er(F o )) 
«int = l^ [|F o - P|(tnean)|]/ J^" 1 



i?32, 155 
9.5203(6) 
7.5439(5) 
592.15(7) 
3.66-38.70:4767 



1/6 
4.712 
774 
137.74 



90(1) K 
GdFe 3 (B0 3 ) 4 
560.04 
Trigonal 

P3i21, 152 
9.5305(3) 
7.5479(2) 
593.73(3) 
3.66-38.55; 54 



0.5 
4.699 
774 
137.37 



ip- and w-scans 
3.66, 38.70 
h : -16 — > 16; k : 
i : -13 -> 13 
0.0563-0.2519 
16.0 
4700 
755 
755 
0.0339 



0.71073 
Graphite 
Bruker SMART APEX, 
CCD area-detector diffractometer 
4096 X 4096/62 X 62 (binned 512) 



ip- and w-scans 
2.47, 38.58 
h : -16 — > 16; k 
I : -13 — > 13 
0.0620-0.2513 
16.0 
13868 
2217 
1956 
0.0338 



c-Reflnement . 



Nu 



of rcflcctk 



Number of refined parameters 
Isotropic secondary-extinction coefficient, g 
Final agreement factors: 

toK(F 2 ) = [^MF,? - Fa ) 2 ]/[™<^r?) 2 ]] 1/2 

Weighting scheme: a, b 

HI = 1/[^ 2 (F 2 ) + (aP) 2 + bP] 

And P = [max(F 2 , 0) + 2F 2 ]/3 

H(F) = ^(||F | - \F c \\)/22 |F |ForF > 4.0cr(F o ) 
Absolute-Structure parameter Flack's x 



0.0371(12) 
0.0374 



-72 \ 2 



r]/(r> 



GooF = S = [2_^[w(F£ - 
n = number of reflections 
p = number of parameters refined 
Residual electron density in final 
Difference Fourier map, e/A3 
Max. (shift/s) final cycle 
Average (shift/er) final cycle 



■ p)] 



1/2 



0.0151 

0.50 

1.075 

-1.0, 0.6(1) 

0.001 
0.000 



0.0359(10) 
0.0452 



0.0312, 0.0 

0.0167 

0.50 

0.855 

-0.5,0.8(1) 

0.003 
0.000 



TABLE II: Fractional atomic coordinates and Anisotropic Parameters at RT 



Atom WyckoiT Schoenflies 


X 


y 


z 


U eq (A) 


Gd 3 a 


D 3 











0.00867(5) 


Fe 9 d 


c 2 


0.21659(5) 


1/3 


1/3 


0.00654(9) 


01 9 e 


c 2 


0.1442(2) 


X 


1/2 


0.0083(3) 


02 9 e 


c 2 


0.4087(3) 


X 


1/2 


0.0132(5) 


03 18 f 


Ci 


0.0254(2) 


0.2125(2) 


0.1824(2) 


0.0095(3) 


Bl 3 b 


D 3 








1/2 ' 


0.0066(6) 


B2 9 e 


c 2 


0.5526(3) 


X 


1/2 


0.0077(4) 



Atom 


Uu 


U22 


^33 


U23 


u 13 


u 12 


Gd 


0.00888(8) 


0.00888(8) 


0.00825(9) 


O.OOOO(-) 


O.OOOO(-) 


0.00444(4) 


Fe 


0.00591(12) 0.00635(14) 0.00750(18) 0.00009(8) 0.00004(4) 0.00317(7) 


Ol 


0.0058(5) 


0.0058(5) 


0.0109(7) 


0.0014(3) 


0.0014(3) 


0.0011(6) 


02 


0.0076(6) 


0.0076(6) 


0.0182(11) 


0.0047(5) 


0.0047(5) 


0.0009(8) 


03 


0.0065(5) 


0.0112(5) 


0.0111(5) 


0.0025(5) 


0.0021(4) 


0.0046(4) 


Bl 


0.0058(9) 


0.0058(9) 


0.0081(14) 


0.0000(0) 


O.OOOO(-) 


0.0029(4) 


B2 


0.0070(7) 


0.0070(7) 


0.0093(8) 


0.0007(4) 


0.0007(4) 


0.0036(11) 



TABLE III: Fractional atomic coordinates and Anisotropic Parameters at T=90 K 



Atom Wyckoff Schocnflies 


X 


y 


z 


U eq (A) 


Gd 


3 a 


c 2 


-0.33342(1) 


X 





0.00406(4) 


Fel 


3 a 


c 2 


0.11536(5) 


X 





0.00360(6) 


Fe2 


6 c 


Ci 


-0.21420(6) 


-0.54975(4) 


0.34154(2) 


0.00366(6) 


Ol 


3 b 


c 2 





-0.07819(15) 


1/6 


0.0072(3) 


02 


6 c 


Ci 


-0.5832(2) 


-0.2709(1) 


0.13774(12) 


0.00692(17) 


03 


6 c 


Ci 


-0.1194(3) 


-0.30445(26) 


-0.17980(18) 


0.0057(3) 


04 


6 c 


Ci 


-0.1467(3) 


-0.36234(16) 


0.18479(18) 


0.0058(2) 


05 


6 c 


Ci 


0.4755(2) 


0.1451(2) 


-0.15980(8) 


0.0057(3) 


06 


3 b 


c 2 


0.1877(2) 





5/6 


0.0049(3) 


07 


6 c 


Ci 


-0.5235(3) 


-0.53811(17) 


-0.18533(18) 


0.0056(2) 


B2a 


6 c 


Ci 


-0.4473(4) 


-0.1201(-) 


0.15617(14) 


0.0053(3) 


B2b 


3 b 


c 2 





-0.2223(3) 


1/6 


0.0049(4) 


Bl 


3 b 


c 2 


0.33204(14) 





5/6 


0.0045(5) 



Atom 

Gd 

Fe 

Fe 

Ol 

02 

03 

04 

05 

06 

07 
B2a 
B2b 

Bl 



Un 
0.00409(6) 
0.00355(10) 
0.00349(9) 
0.0081(6 
0.0058(3 
0.0062(4 
0.0054(4 
0.0051(4 
0.0044(4 
0.0047(4 
0.0050(6 
0.0018(8 
0.0042(8 



U22 
0.00409(7) 
0.00355(10) 
0.00356(11) 
0.0054(4) 
0.0062(3) 
0.0040(5) 
0.0045(4) 
0.0050(4) 
0.0058(6) 
0.0063(4) 
0.0067(7) 
0.0061(6) 
0.0048(9) 



U33 
0.00408(7) 
0.00355(10) 
0.00378(9) 
0.0089(4) 
0.081(3 
0.0065(4) 
0.0067(4) 
0.0063(5) 
0.0062(5) 
0.0064(4) 
0.0055(5) 
0.0055(7) 
0.0046(11) 



U23 
0.00010(2) 
0.00006(4) 
0.00013(7) 
0.00121(19) 
0.0012(3) 
0.0006(3) 
0.0007(3) 
0.0007(3) 
0.0019(5) 
0.0008(3) 
0.0004(5) 
0.0004(3) 
0.0009(5) 



Ua 
0.00009(1) 
0.00006(4) 
0.00003(5) 
0.0024(4) 
0.0014(3) 
0.0007(3) 
0.0003(3) 
0.0003(3) 
0.0010(2) 
0.0001(3) 
0.0001(4) 
0.0008(6) 
0.0004(2) 



U12 
0.00210(3) 
0.00166(10) 
0.00166(6) 
0.0040(3) 
0.0025(3) 
0.0023(4) 
0.0019(4) 
0.0019(5) 
0.0029(3) 
0.0033(4) 
0.0040(4) 
0.0009(4) 
0.0024(4) 



nearly perpendicular to the C3 axis and the Fe atoms are arranged in helicoidal chains parallel to this axis. Different 
chains arc connected by GdC>6 and BO3 groups, where each individual BO3 and GdC-6 group connects three chains. 
The distance between Fe atoms in the same chain (3.1669(4) A) is shorter than the Fe-Fe distances for two nearest 
chains, that varies along the chain, 4.8308(5) A, being the shortest one. The main exchange interaction between Fe 3+ 
is therefore of quasi- ID nature. 



TABLE IV: Symmetry position of BO3 groups, angle between the group and the c-axis and their flatness of them The flatness 
is expressed as distance of the B atoms from the plane defined by the three oxygen ligands. 

Room Temperature Low Temperature 

Symm. Angle [degr.] Flatness Symm. Angle [degr.] Flatness 
Bl D 3 90 Bl C 2 87.52(4) 0.00002(4) 

B2 C 2 84.37(11) B2a Ci 81.89 0.0055(11) 

B2b C 2 83.55 0.0000(2) 

Upon lowering the temperature the GdFe3(B03)4 crystal reduces the symmetry from R32 to P3i21, in the trigonal 
system. Fig. [SJshows the coordination polyhedra (GdO@, FeOe, and BO3) for two different structures i?32 and P3i21. 
At RT the BO3 groups occupy two inequivalent positions, Bl (D3) and B2 (C2). At 90 K the site symmetry of the Bl 
atoms is reduced to C2, whereas the site symmetry of the B2 atoms differentiates into a B2b (C2) and B2a (Ci) (see 
also Table Hvjl . In the LT-phase, the angle between BO3 groups and c-axis is changed. Moreover, in the LT-phase 
the BO3 groups in Ci position are considerably distorted and no longer flat (see Tab lIV|) . 

The LT structure is in agreement with the Raman datai^. The group theoretical vibrational analysis, based on the 
90 K structure, shows that, due to the lowering of the symmetry, new librational modes of the BO3 become Raman 
active (R x and R y ). The observed structural change is compatible with an appearance of the librational R y mode 
of the BO3 group upon approaching T s from above, and a subsequent hardening of this mode in the LT phase (See 
FiglHJl. Also the anisotropic displacement parameters (Table ITTfl L/33 of the oxygen atoms suggest that the borate 
groups are relatively free to oscillate around the y-axis. 

Concerning the magnetic structure, the main peculiarity of both the RT and the 90 K structures is the existence of 




FIG. 2: Coordination polyhedra for the RT (7?32-(a)) and the LT (P3i21-(b)) structures of GdFe 3 (B0 3 )4- At RT the Fe atoms 
with C2-symmetry are coordinated by three different types of oxygen atoms. AT 90 K there are two inequivalent positions 
for the Fe atoms: Fel (the same symmetry as RT) and Fe2, which is in a general position (surrounded by five different types 
of oxygen). At RT the Gd atoms are surrounded by six oxygens of one type (D3-symmetry). At 90 K they are coordinated 
by three different types of oxygen atoms (C2-symmetry). For the B atoms there are two kinds of coordination at RT. Bl (in 
D3-symmetry) is surrounded by three oxygen atoms of the same type, the BlOa-group is thus an equilateral triangle. The 
B2 atoms are surrounded by two types of O, B2O3 is an isosceles triangle (C2-symmetry). At 90 K there are three kinds of 
coordinations for the B atoms: B2b03 and Bf03 are isosceles triangles (C2 symmetry) and B2a03 is general triangle. 

magnetically quasi-lD helicoidal iron chains (see Fig. The intra-chain exchange interaction between the Fe ions 
is expected to be dominated by Fe-Fe direct exchange and Fe-O-Fe superexchange, depending respectively on Fe-Fe 
distance and two Fe-O-Fe angles (FigHJ. Nevertheless there are some significant differences. 

At RT all the Fe atoms are in equivalent positions (C2). All Fe-chains are equivalent as are the Fe-O-Fe angles 
(102.40°(12) and 103.65°(8)) and Fe-Fe distances (3.1669(4) A). Therefore, the exchange interactions between neigh- 
boring iron ions within a chain are also equivalent. At LT (Fig. QJ, as shown in Table ITVl the BO3 groups form 
one general triangle (gray Ci position) and two isosceles triangles (white and dark, C2 position). At RT "gray" and 
"white" groups become also equivalent, with C2 symmetry, and the "dark" one is a regular triangle with the D3 
symmetry. The reduction of the symmetry of the borate groups changes the surrounding of the Fe atoms (Fig. |5J), 
yielding two inequivalent positions (C2 and Ci). Therefore, the Fe-Fe distances are different for the two chains: one 
is stretched (3.1828(4) A) and in the other one is compressed (3.1554(4) A). The angles Fe-O-Fe for the first chain are 
101.24°(5) and 103.71°(9), while those for the second are 102.46°(6) and 103.91°(8). Therefore also the intra-chain 
exchange interaction is different for the two chains. 

The GdOe prism connects three chains, one containing Fcl-atoms with C2 symmetry and two with Fe2-atoms in 
a general position. There is only one inequivalent position for the Gd-ions in both the RT and the 90 K structures. 
However, the site symmetry of the Gd-ions changes from D3 at RT to C2 at 90 K. This confirms the interpretation of 
the authors of Refi, who interpreted low-temperature infrared spectra of Ndo.oiGd .9gFe3(B03)4 in terms of Kramers 
doublets of Nd, assuming only one structural position for the Nd ions. 

The detailed structure of the Fe-chains and their interconnection are shown in Fig. ^ The intra-chain interactions 
between Fe-atoms go through two Fe-O-Fe superexchange pathways. The magnetic inter-chain interaction arises 
from the Fe-O-O-Fe, and possibly the Fe-O-Gd-O-Fe superexchange path. The role of the first superexchange path 
is important because the substitution of Gd for non-magnetic Y does not lead to a disappearance of 3D magnetic 



50 

E 
o 

I 40 

'</> 
o 

Q- 

CD 

5 30 



50 100 150 

Temperature (K) 

FIG. 3: Temperature dependence of Raman frequency of the low energy mode appearing at the structural phase transition T s . 
Insert shows the low-frequency part of Raman spectrum at two close temperatures before and after phase transition. 




FIG. 4: Two Fe-chains and the most important exchange paths: the intra-chain exchange is via Fe-Fe direct exchange or 
Fe-O-Fe super-exchange, while the inter-chain exchange is through Fe-O-Gd-O-Fe and Fe-O-O-Fe. 

ordering. On the contrary, the Neel temperature for YFe3(B03)4, Tni=40 K, is larger than the one found for Gd- 
compound (T N i=37 K), and for NdFe 3 (B0 3 ) 4 (T N i=3Q K)2J±. Considering the different RE 3+ ionic radii (0.983 A 
for Nd, 0.938 A for Gd, and 0.900 A for Y) it is clear that the Neel temperature of RFe3 (1303)4 depends strongly 
on the ionic radii: a smaller ionic radius results in a higher Tni- Moreover, Tni does not seem to be affected by 
the spin of the RE. In this sense it is clear that the main superexchange path of the interaction between different 
chains is the Fe-O-O-Fe path, and that a small distortion of this path changes substantially the magnetic properties 
of the system. It is therefore clear that the existence of two nonequivalent Fe chains with different Fe-Fe distances 
and Fe-O-Fe angles could lead to substantially different intra-chain exchange constant for the two chains^, therefore 
it should be taken into account for the interpretation of the magnetic properties. 

IV. CONCLUSIONS 

In summary, we have determined the crystal structure of GdFe 3 (603)4 at RT and 90 K. At RT GdFe3 (603)4 
exhibits R32 structure, in agreement with Refsi^. Below the structural phase transition (T s = 156 K) the structure 




has the P3\21 space group. The main difference of LT-structure compared to the RT one is the lowering of the 
symmetry and the tilt of the BO3 groups. This confirms the interpretation of the Raman spectra of this phase 
transition. The main conclusion resulting from the LT structure determination is the presence of two inequivalent 
positions for the Fe atoms giving rise to two different iron hclicoidal chains. 
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